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ABSTRACT 

We present constraints on cosmological parameters based on a sample of Sunyaev-Zel'dovich-selected 
galaxy clusters detected in a millimeter-wave survey by the Atacama Cosmology Telescope. The cluster sam- 
ple used in this analysis consists of 9 optically-confirmed high-mass clusters comprising the high-significance 
end of the total cluster sample identified in 455 square degrees of sky surveyed during 2008 at 148 GHz. We 
focus on the most massive systems to reduce the degeneracy between unknown cluster astrophysics and cos- 
mology derived from SZ surveys. We describe the scaling relation between cluster mass and SZ signal with 
a 4-parameter fit. Marginalizing over the values of the parameters in this fit with conservative priors gives 
er 8 = 0.851 ±0.115 and w = -1.14 ±0.35 for a spatially-flat wCDM cosmological model with WMAP 7-year 
priors on cosmological parameters. This gives a modest improvement in statistical uncertainty over WMAP 
7-year constraints alone. Fixing the scaling relation between cluster mass and SZ signal to a fiducial relation 
obtained from numerical simulations and calibrated by X-ray observations, we find a% = 0.821 ±0.044 and 
w = -1 .05 ± 0.20. These results are consistent with constraints from WMAP 7 plus baryon acoustic oscillations 
plus type la supernoava which give erg = 0.802 ±0.038 and w = -0.98 ±0.053. A stacking analysis of the 
clusters in this sample compared to clusters simulated assuming the fiducial model also shows good agreement. 
These results suggest that, given the sample of clusters used here, both the astrophysics of massive clusters and 
the cosmological parameters derived from them are broadly consistent with current models. 
Subject headings: cosmic microwave background - galaxies: clusters: general - cosmology: observations 
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1. INTRODUCTION 

Ever-improving observations suggest a concordant picture 
of our Universe. In this picture, generally called ACDM, 
"dark energy," the component responsible for the Universe's 
accelerated expansion, is believed to be the energy of the vac- 
uum with a constant equation of state parameter, w, equal 
to -1 (e.g.JRiess et alJl2009l:lBrown et al]|2009tlHicken et al] 
2009t [Kessler et alj|2009t iPercival et alj|2010t iKomatsu et alj 



2010). ACDM has been measured via probes of the Uni- 
verse's expansion rate such as type la supernovae, the pri- 
mary cosmic microwave background, and baryon acoustic 
oscillations. However, ACDM also makes concrete predic- 
tions about the Universe's growth of structure. This growth 
rate describes how quickly dark matter halos form and evolve 
over cosmic time. A deviation from this predicted growth 
rate, particularly on li near sca l es, wo u ld signal a breakdown 
of ACDM (see e.g.. funded (120051): iBertschinger & Zukir] 
d2008b: ISilvestri & TroddeTd2009h : llain & Khourvl d2010l) : 
Shapi ro et alj d2010h and references therein). 

A handful of techniques have been available for measur- 
ing the growth of structure in the Universe. These largely 
consist of observing the weak and strong lensing of back- 
groun d sources by intervening matter (e.g., Schrabbac k et alj 
2010), measuring distortions in redshift space with spectro - 
scopic surveys of galaxies (e.g., Simp son & Peacockll2010l) . 
and quantifying the abundance of galaxy clusters as a func- 
tion of mass and redshift (e.g., Eahcall & Fan 1998)). This 
latter technique is one of the oldest an d has been ma turing 
with the adv ent of large area X-r ay (e. g.. iTruemperl 1 990h and 
optical (e.g jKoester et a l. 2007) surveys. 

Millimeter-wave surveys now possess the resolution and 
sensitivity to detect galaxy clusters. Detecting galaxy clus- 
ters via the Sunyaev-Zel'dovich (SZ) effect in large area 
millimeter-wave maps, as have become a vailable through the 
Atacama Cosmology Telescope ( ACT) (|Swetz et aT] 120101) 
and the South Pole Telescope (SPT) dCarlstrom et al.l2009l) . is 
a poten tially powerful method. Cluster selection using the SZ 
effect (IZei'dovich & Sunyaevl Il969t ISunyaev & Zel'dovichl 
1 19701 119721) is the technique whose selection function is least 
dependent on cluster redshift. This allows for a complete pic- 
ture of the evolution of clusters from their first formation to 
the present. 

Here we probe structure growth with a measurement of the 
abundance of massive galaxy clusters from observations made 
by the ACT project in 2008. We focus on the most massive 
SZ-selected clusters as this is the regime where high signal- 
to-noise measurements exist and we can best understand the 
cluster astrophysics. We also note that the clusters consid- 
ered in this work are rare and represent the tail of the mass 
distribution, which is sensitive to the background cosmol- 
ogy. This analysis uses the number of massive galaxy clus- 
ters to constrain, in particular, the normalization of the matter 
power spectrum, cs, and the dark energy equation-of state- 
parameter, w. 

This paper is structured as follows: Section|2]describes the 
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SZ effect and the ACT SZ cluster survey. Section[3]describes 
the 2008 ACT high-significance cluster sample. In Section 
|4j we present our results, and in Section [3] we discuss their 
implications and conclude. 

2. BACKGROUND 

2.1. The Thermal Sunyaev-Zel'dovich Effect 

The thermal SZ effect arises when primary cosmic mi- 
crowave background photons, on their path from the last 
scattering surface, encounter an intervening galaxy cluster. 
The hot ionized gas within the cluster inverse Compton 
scatters about 1% of the CMB photons, boosting their en- 
ergy and altering the intensity of the microwave background 
as a function of frequency at the location of the cluster 
dSunvaev & Zel'dovicrJ[T97all972h . To first order, the effec- 
tive temperature shift (which is proportional to the intensity 
shift), at a frequency v, from the thermal SZ effect is given by 



Ar tl \ kB(JT I T V 
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where n e and T e are the number density and temperature of 
the electron distribution of the cluster gas, dl is the line-of- 
sight path length through the cluster, <tt is the Thomson cross- 
section, kg is the Boltzmann constant, and 



f(x) = xcoth(x/2)-4, x = hv/(k B TcMB)- 



(2) 



Here y is the usual Compton y-parameter. No te that the full 
SZ eff ect contains relativistic corrections as in iNozawa et al.l 
( 1998), which we include in our simulations. For our sample, 
these corrections are 5 to 10%. We take Eq. [T] describing 
the first-order thermal SZ effect, as the definition of y we use 
throughout this work, and treat the relati vistic corrections as 
an additional source of noise (see Section [3~3l ). 

At frequencies below 218 GHz, where the signal is null, 
AT is negative, and the cluster appears as a cold spot in CMB 
maps. Above the null, AT is positive, and the cluster ap- 
pears as a hot spot. Eq.[T]is also redshift independent, and the 
amplitude of the intensity shift is to first order proportional 
only to the thermal pressure of the cluster. This makes the 
SZ effect especially powerful for two reasons: the microwave 
background can trace all the clusters of a given thermal pres- 
sure that have formed between the last scattering surface and 
today in a redshift independent way, and the amplitude of this 
effect, being proportional to the thermal pressure, is closely 
related to the cluster mass. 

2.2. The ACT Sunyaev-Zel 'dovich Cluster Survey 

The Atacama Cosmology Telescope (ACT) is a 6-meter off- 
axis telescop e designed for arcminute-scale m illimeter-wave 
observations dSwetz et al.ll2010l : lHincks et al.ll2009l) . It is lo- 
cated on Cerro Toco in the Atacama Desert of Chile. One 
goal of this instrument is to measure the evolution of struc- 
ture in the Universe via the SZ effect. In the 2008 observ- 
ing season ACT surveyed 455 square degrees of sky in the 
southern hemisphere at 148 GHz. In this surv ey, galaxy clus- 
ters we re detected from their SZ signal (see Marriage et al. 
(2010a) for details). A sample of 23 SZ-selected clusters 
was optically confirmed using multi-band optical imaging 
on 4-meter telescopes durin g the 2009B observing season 
(see iMenanteau etafl d201 0a) for details). Some of the low- 
redshift systems in this sample are previously known clus- 
ters for which spectroscopic redshifts are available. However, 
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roughly half are newly detected systems, and photometric red- 
shift estimates have been obtained from optical imaging. Here 
we make use of the subsample of these clusters with high- 
significance SZ dete ctions (signal-to-noise ratio > 5, as de- 
fined in Section [3T2b to obtain cosmological parameter con- 
straints. 

3. THE 2008 ACT HIGH-SIGNIFICANCE CLUSTER 
CATALOG 

3.1. CM B Data 

Here we give a brief overview of the survey observations 
and the reduction of the raw data to maps. For a more com- 
plete introduction to the ACT instrument, obse rvations, and 
data re duc tion pipeline , we r efer the reader to Fowler et al. 
(2010) and lSwetz etail (120101) . The 2008 observations in the 
southern hemisphere were carried out between mid-August 
and late December over a 9° wide ACT strip centered on 
a declination of -53° degrees and extending from approxi- 
mately 19 h through A to 7 /! 30 in right ascension. The 455 
square degrees used for this analysis consists of a 7°-wide 
strip centered at a declination of -52°30' and running from 
right ascension 00 ft 12 m to 7 ft 10 m . The resolution of the ACT 
instrument is about 1.4' at 148 GHz. Typical noise levels in 
the map are 30 /iK per square arcminute, rising to 50 /iK to- 
ward the map boundaries. Seven of the nine clusters consid- 
ered in this work fall in the central region of the map with 
lower noise levels. 

Rising and setting scans cross-link each point on the sky 
with adjacent points such that the data contain the information 
necessary to make a map recovering brightness fluctuations 
over a wide range of angular scales. In addition to survey ob- 
servations, ACT also executed regular observations of Uranus 
and Saturn during 2008 to provide beam profiles, pointing, 
and temperatur e calibration . Ana lysis of the beam profiles 
is discussed in Hinc ks et al.l (120091) . Absolute pointing is de- 
termined by comparing the positions of ACT-observed radio 
sources with the pos itions of these sam e sources detected in 
the AT20G survey dMurphv et al.ll2Q10t) . The final tempera- 
ture calibration at 148 GHz is based on a recent analysis cross- 
correlating ACT and WMAP maps and is determined with an 
uncertainty of 2% (Hajian et al. 2010). The small residual 
calibration uncertainty translates into a small systematic un- 
certainty in y values of observed clusters, which is negligible 
compared to other uncertainties discussed in this analysis. 

To make maps, an iterative preconditioned conjugate gradi- 
ent solver is used to recover the maximum likelihood maps. 
This algorithm solves simultaneously for the millimeter sky 
as well as correlated noise (e.g., a common mode from at- 
mospheric emission). The map projection used is cylindrical 
equal area with a standard latitude of -53°30' and square pix- 
els, 0.5' on a side. 

3.2. Cluster Detection Method 

In order to detect clusters in single-frequency millimeter- 
wave maps we construct a filter that is similar in morphology 
to the clusters we are trying to detect. We adopt a matched 
filter of the form 

^ V [{lit) 2 J P(k>) J P(k) 

followi ng lHaehnelt & Tegmarkl (11996b . iHerranz et al.l 
d2002allrl . and iMelin et alj (l2006l) . Here r(k) is the beam 



convolved cluster signal in Fourier space, and P(k) is the 
power spectrum of the noise, both astrophysical and instru- 
mental. The astrophysical noise sources for cluster detection 
include the primary CMB lensed by intervening structure, 
radio galaxies, dusty star-forming galaxies, Galactic dust, and 
the SZ background from unresolved clusters, groups, and the 
intergalactic medium. Since the power from the SZ signal is 
subdominant to the se astrophys i cal sources (as evidenced by 
Lueker etalj (120101) . lHall et all (120 10). Fo wler et al l d2010l) . 
Das et al] (120101) . and iDunkley et al.l (120101) ). we can to a 
good approximation model the power spectrum of the total 
noise as the power spectrum of the data itself. In Eq. [3] the 
quantity in square brackets serves as a normalization factor 
to ensure an unbiased estimate of the cluster signal. When 
multi-frequency maps are available, this filter can be modified 
to incorporate the known spectral signature of the SZ signal. 

The template shape that we choose to match the cluster 
morphology is given by a two-dimensional Gaussian profile, 
which in Fourier space has the forrrQ 

AT(l)=Aexp[-9 2 (l + 1)1/2]. (4) 

Here 9 = FWHM/VShd where FWHM is the fuU width at 
half maximum, and A is a normaliza tion factor that will be de- 
rived from simulations (see Section [373l . We choose FWHM 
to be 2' as this is a typical cluster size in our maps. The analy- 
sis presented here of the cosmological parameters is nearly in- 
dependent of the particular profile chosen for the cluster tem- 
plate, as long as the template is smooth and well-matched to 
the cluster angular size. 

Before filtering our map to find clusters, we multiply the 
map, pixel-wise, by the square root of the number of obser- 
vations per pixel normalized by the observations per pixel in 
the deepest part of the map in order to establish uniform noise 
properties. We then detect point sources (radio and infrared 
galaxies) by a matched filter with the ACT beam as the tem- 
plate. Selecting all point sources with a signal-to-noise ratio 
greater than 4.0 in this filtered map, we mask them by re- 
placing all on-source pixels with signal-to-noise ratio greater 
than 4.0 with the average of the brightness in an annulus 4' 
away from the source center. We do this to avoid false de- 
tections due to the filt er ringing around bright sources. See 
iMarriage et al] d2010bl) for details regarding point source de- 
tection. 

After masking out the brightest point sources, we filter the 
map to find clusters. Clusters are then detected within this 
filtered map with a sim ple peak detection algo rithm along 
the lines of SExtractor (Ber lin & Arnoutsl [l996). An SZ y 
value for the brightest 0.5' pixel is measured for each clus- 
ter using this filtered map. This definition of y is different 
from the integrated Y, which is specifically the Compton-y 
parameter integrated over the face of the cluster and given by 
Y = f yd ft. The integration for this Y value is performed over 
a radius tied to the size of the cluster, and a Y defined this 
way would be a preferable qu antity to use, having lower scat- 
ter with mass in theory (e.g.. Ida Silva et al]|2004t iMotl et al] 
l2005llNagail2006l:lReid & Spergell2006HBonaldi et alj2007l) . 
However, given single-frequency millimeter-wave maps, the 
size of each cluster cannot always be robustly determined. An 
alternative quantity to measure is a "central y value," gener- 
ally referred to as yo> which essentially describes the normal- 
ization of the specific template shape used to find the cluster. 

1 We use the flat space approximation, 1 = Ink. 
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TABLE 1 

ACT Cluster Catalog for High-Significance Clusters from the 2008 Observing 

Season 



ACT Descriptor 



R.A. 



decl. 



.yTcMB(^K) 



Redshift 



Other Name 



ACT-CL J0645- 


-5413 


06:45:30 


-54:13:39 


340 ± 60 


0.1 67 a 


Abell 3404 


ACT-CL J0638- 


-5358 


06:38:46 


-53:58:45 


540 ± 60 


0.222 a 


Abell S0592 


ACT-CL J0658- 


-5557 


06:58:30 


-55:57:04 


560 ± 60 


0.296 b 


lES0657-558(Bullet) 


ACT-CL J0245- 


-5302 


02:45:33 


-53:02:04 


475 ± 60 


0.300 c 


Abell S0295 


ACT-CL J0330- 


-5227 


03:30:54 


-52:28:04 


380 ± 60 


0.440 d 


Abell 3128(NE) 


ACT-CL J0438- 


-5419 


04:38:19 


-54:19:05 


420 ± 60 


0.54±0.05 c 


New 


ACT-CL J0616- 


-5227 


06:16:36 


-52:27:35 


360 ± 60 


0.71±0.10 c 


New 


ACT-CL J0102- 


-4915 


01:02:53 


-49:15:19 


490 ± 60 


0.75±0.04 c 


New 


ACT-CL J0546- 


-5345 


05:46:37 


-53:45:32 


310 ±60 


1.066 f 


SPT-CL 0547-5345 



^ /iK given for the brightest 0.5' pixel of each cluster 

a spectroscopic-z from de Grandi et al. (1999) 

k spectroscopic-z from Tucker et al. (1998) 

c spectroscopic-z from Edge et al. 1 1994) 

^ spectroscopic-z from Werner et al. ( 2007) 

e photometric-z from Menanteau et al. 1 20 10a) 

^ spectroscopic-z from Infante et al. ( 2010); Brodwin et al. 1 2010) 



This quantity is not ideal for a cosmological analysis as it is 
intimately tied to the profile shape whereby the yo value is 
determined. These values also exhibit a larger s catter with 
cluste r mass than an integrated y quantity (e.g., iMotl et all 
120051) . For all the clusters considered in this analysis, we fix 
an aperture size, given by our pixel size of 0.5', and mea- 
sure Compton-y values within this fixed aperture. We do not 
consider larger aperture sizes here because we wish to limit 
template shape and redshift dependence. 

We find that selecting clusters with a yTcMB valu^l greater 
than 300 /iK corresponds to a subsample of clusters with 
a signal-to-noise ratio greater than 5. Here signal-to-noise 
ratio is defined as the signal of the brightest cluster pixel 
in the filtered map divided by the square root of the noise 
variance in the filtered map. This subsample corresponds 
to the subsample o f cluste rs with signal-to-noise ratio > 5.9 
in iMarriage et alj (l2010al) PI The y values derived with the 
m ethod used in th i s paper are not directly comparable to those 
in IMar riage et al. (2010a), in which an optimal filter sized is 
searched for. However, the methods are independently com- 
pared to simulations. While we detect clus ters down to a 
signal-t o-noise ratio of about 3 as defined in Marriage et al. 
(2010a), we use only this higher-significance subsample in 
this work. This subsample is given in Table [1] 

3.3. Simulations and SZ Signal Recovery 

To determine the expected scatter in our recovered y7cMB 
values, we perform the same detection procedure discussed 
above on simulated maps. H ereafter, the s imulat ions we re- 
fer to are those discussed in Sehg al et al.1 (120101) . which in- 

2 Note that y is a dimensionless parameter. We multiply it by 7cmb = 
2.726 X 10 6 /^K to give an indication of the expected temperature decrements. 
For the frequency dependence, f(x) —1 in Eq.[T]at 148 GHz. 

3 In Marriage et al. 1 2010a), a different detection method is used that varies 
the angular scale of the filter to match clusters of different sizes, and assigns 
a signal-to-noise ratio based on the scale that gives the highest value. The one 
cluster that has a signal-to-nois e ratio > 5.9 in that wo rk that is not included 
here is ACT-CL J0235-5121. In Marriage efaTl KOlOal) . it was found to have 
a high signal-to-noise ratio using a template scale of 4.0' . Although there is 
no doubt that this is a massive cluster (Menanteau et al. 2010a), its relatively 
high redshift, z = 0.43 ± 0.07 argues for a compact size, suggesting that CMB 
contamination could b e boosting the clusters 's signal-to-noise ratio on a 4.0' 
scale, as discussed in Marriage et al. 1 2010a). This cluster is not found with 
signal-to-noise ratio > 5 using the 2' FWHM Gaussian template described 
above. 



elude the SZ signal, lensed primary cosmic microwave back- 
ground, Galactic dust, and radio and infrared sources cor- 
related with SZ clusters as suggested by observations. The 
large-scale structure in thi s simulation was carried out using 
a tree-particle-mesh code (Bode et al. 2000; Bode & Ostrikerl 
2003), with a simulation volume of 1000 /i _I Mpc on a side 
containing 1024 3 particles. The cosmo logy adopted is con- 
sistent with the WMAP 5-year results (iKomatsu et al.ll2009l) 
though the details of the cluster properties are relatively in- 
sensitive to the background cosmology. The mass distribution 
covering one octant of the full sky was saved, and halos with a 
friends-of-friends mass above 1 x 10 I3 M Q and with a redshift 
below z = 3 are identified. The thermal SZ signal is derived by 
adding to the N-body halos a gas prescription that assumes a 
polytropic equation of state and hydrostatic equili brium. Thi s 
model, which is described in more detail in lBode et al.l (12009). 
adjusts four free parameters (star-formation rate, nonthermal 
pressure support, dynamical energy transfer, and feedback 
from active galactic nuclei) which are calibrated against X- 
ra y gas fractions as a fu nction of temperature from the sample 
of ISun et al.l (120091) and lVikhlinin et al.l (120061) . The pressure 
profiles of the massive, low-redshift cl usters in this simula- 
tion agree well with the best-fit profile of lArnaud et al.1 (120091) 
based on X-ray obse rvations of high-mass, low- reds hift sys- 
tems dTrac etal.ll2010l) . We will see in Section 03] that the 
stacked SZ signal of the clusters in Table Q] is also consistent 
with the stacked thermal SZ signal of the massive clusters in 
this simulation. The kinetic SZ in this simulation is calculated 
from the line-of-sight momentum of the particles. We also 
in clude the relat i vistic corrections to the SZ signal as given 
in lNozawa et al.l dl998l) . We convolve these simulations with 
the ACT beam and run them thro ugh the same map-making 
process discussed in Section |3~T1 including simulated atmo- 
spheric emission and realistic instrumental noise. 

From these simulations, we cut out six different patches of 
455 square degrees to mimic the sky coverage in this analysis. 
These six sky patches give us about 40 clusters that would cor- 
respond to the high-significance cluster sample given in Table 
[U Using these simulations, we apply the same cluster detec- 
tion procedure as discussed in Section l3~2l and recover y7cMB 
values for the detected clusters. These recovered yTcuB val- 
ues are compared to the true yTcMB values taken from the 
first-order thermal SZ maps alone, prior to any instrumental 
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FIG. 1 . — The true y value (as denned in Eq.[T} for the brightest 0.5' pixel 
of each cluster versus the recove red y value from simulations using the detec- 
tion method outlined in Section |X2| The dimensionless y values have been 
multiplied by 7cmb = 2.726 X 10 6 to give an indication of the expected 
temperature decrements at 148 GHz. The root-mean-square scatter shown 
here is 60 /iK. 

or atmospheric modifications. The comparison between these 
true and recovered v7cmb values is shown in Figure Q] We 
set the normalization factor, A, in Eq. [4] such that the mean 
bias between true and recovered y7cMB values is zero. The 
root-mean-square scatter in the recovered yTcMB values is 60 
/iK. The scatter here is dominated by the instrumental and 
atmospheric noise sources in the map and not by the tech- 
nique itself. This same normalization is used for filtering the 
data map, and the same scatter is assumed. Figure [TJ also 
shows that any bias due to boosting a cluster with an intrinsic 
yTcMB < 300 /xK to a va l ue a bove 300 [iK is far below the 
scatter. Vanderl inde et al.l ( 120101) found that this signal boost- 
ing effect is at most 4% for clusters detected with a signal-to- 
noise ratio greater than 5. 

3.4. Optical Identification and Cluster Redshifts 

The sample of SZ-selected clusters obtained from the data 
via the method above was followed up with optical obser- 
vations to verify the millimeter-wave cluster identifications 
and determine cluster redshifts. Here we provide a summary 
of the observing strategy and redshift deter minations of our 
cluster sample, and we refer the reader to Me nanteau et al.l 
(2010a) for a detailed description. 

SZ cluster candidates were observed during the 2009B ob- 
serving season with optical imaging on the 4-meter SOAR and 
NTT telescopes to search for a brightest cluster galaxy and an 
accompanying red sequence of cluster members. While some 
of the clusters in the SZ-selected sample correspond to previ- 
ously known systems at low redshift (z <0.3), some represent 
new systems, previously undetected at other wavelengths. 

Photometric redshifts and their probability distributions, 
p(z), were computed for each object from thei r dust-correcte d 
gri isophotal magnitude using the BPZ code (Bemtez 2000). 
Six of the clusters in our sample have spectroscopic redshift 
information available (see Table Q] for references), as they 
were previously known systems. For the remaining three sys- 
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FIG. 2. — The redshift distribution of the 9 high-significance c lusters listed 
in Table [T] compared to expectations from the mass function of Tin ker et alj 
12008) assuming all clusters above a given mass threshold have been detected. 
Here we assume the total area observed is 455 square degrees with uniform 
coverage, and the WM AP7+BAO+SN best-fit cosmology for a wCDM model 
I Komatsu et al. 2010) for the mass function. This figure suggests an effective 
mass threshold of our sample of ss 10.4 X 1O 14 M0 (dashed line). For illus- 
trative purposes we show via dotted lines the expected redshift distributions 
for mass thresholds larger/smaller by 0.6 X 10 i4 Mq. 

tems we provide photometric redshifts based on the NTT and 
SOAR imaging. Table[T]gives the mean photometric redshift 
for these clusters obtained by iteratively selecting galaxies 
photometrically classified as E or E/SOs within a projected 
radius of 500 kpc and redshift interval |Az| = 0.05. This 
was done to obtain a local color-magnitude relation for each 
color combination using a 3a median sigma-clipping algo- 
rithm. The uncertainties on the photometric redshifts of the 
three new clusters in Table Q] come from the weighted rms 
of the individual galaxies cho sen as members. Details of the 
photo-z algorithm are given in Menante au et alj (|2010b). 

4. RESULTS 

4.1. SZ Selection Function 

We determine the selection function of our cluster subsam- 
ple through both optical observations and simulations. To 
investigate the effective mass threshold of our cluster sam- 
ple, we plot cumulative clusters as a function of redshift 
and compare that to expectations from the mass function of 
iTinker et a l. (2008) assuming the best-fit cosmology from 
WMA P7+BAO+SN with a wCDM model dKomatsu et alj 
2010). We find that this sample is consistent with a mean 
masf| threshold of 10.4 x 10 14 M Q for M 2 qo as shown by the 
dashed line in Figure [2] For illustrative purposes, we show 
by dotted lines the expected redshift distributions for mass 
thresholds larger and smaller by 0.6 x 10 14 M Q , which enclose 
our subsample distribution for z > 0.25. This figure also as- 
sumes uniform coverage over an area of 455 square degrees. 
Note that Figure [2] is presented to give a qualitative under- 

4 Note that throughout this text cluster masses are defined in terms of M200, 
which is the mass within Rjqo, the radius within which the mean cluster den- 
sity is 200 times the average density at the cluster redshift. 
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Fig. 3. — Completeness versus true v7cmb from simulations. For an ob- 
served yTcMB threshold of 300 fiK we expect our sample to be about 85% 
complete for true \>7cmb values above 300 /xK. 

standing of our sample, but is not directly used in the cosmo- 
logical analysis. 

We also use simulations to characterize our sample . Th e 
same cluster detection procedure discussed in Section [3~2l is 
applied to th e simulated millimeter-wave maps discussed in 
Section 13.31 including instrumental and atmospheric noise 
sources. We find for an observed yTcMB threshold of 300 /iK, 
we expect from simulations to be about 85% complete above 
a true yTcMB value greater than 300 fiK as shown in Figure 
[3] This completeness is calculated as the total number of ob- 
served clusters with a recovered y value above the threshold 
versus the total number of expected clusters with a true y value 
above the threshold. This takes into account the scatter of ob- 
served y values across the threshold. We take 300 /iK as our 
threshold yTcMB value and only consider clusters with mea- 
sured yTcMB values larger than this. We expect about 90% 
purity for cluster detections with an observed yTcMB value 
greater than 300 /iK. Figure |4]illustrates how this false detec- 
tion rate is expected to vary as a function of observed yTcMB 
threshold. From the optical observations discussed above, we 
find that all clusters identified in the millimeter-wave maps 
with yTcMB values greater than 300 /iK were verified as clus- 
ters in the optical. Thus this sample is 100 % pure with no 
false detections (see Men anteau et al.ld2010al) ). Note that Fig- 
ures|3]and|4]are representative, coming from a relatively small 
sample of si mulated high-significance clusters as discussed in 
Section 13.31 They are presented to give a qualitative under- 
standing of the sample and to understand above what yTcMB 
threshold our sample is roughly complete. Beyond this, they 
do not enter in the analysis of cosmological parameters. 

4.2. Scaling Relation Between SZ Signal and Mass 

We assume for the relation between the SZ signal and mass 
the general parameterized forrrQ 

M£)'(£) C ' °> 

s Note that the f orm of this relation is analogous to that in 
IVanderlindeeTaTI fSCTTOh . 
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Fig. 4. — Purity versus recovered }'7cmb from simulations. For an observed 
y7cMB threshold of 300 fiK we expect our sample to be about 90% pure. Our 
actual sample is 100 % pure, as each ca ndidate cluster has been confirmed by 
optical observations i Menanteau et al. 2010a). 

where M = 5 X 10 U M Q hr l and z = 0.5. We also assume 
that cluster y values are randomly distributed around this 
relation with a lognormal scatter S. By lognormal scatter, 
we mean that the scatter in the relation ln(y true ) = ln(A) + 
Bln(M true /M ) + Cln((l +z)/(l +z )) has a Gaussian distribu- 
tion with a mean of zero and a standard deviation equal to S. 
This scaling relation relates the observable quantity for each 
cluster, SZ signal, to the quantity most directly relevant for 
cosmology, the cluster mass. This form is chosen as the SZ 
signal is expected to have a power-law dependence on mass. 
Some dependence of our recovered SZ signal on redshift is al- 
lowed with an additional parameter C. Note that the B and C 
parameters are independent of Mo and zo, whose values only 
affect A. We also note that the uncertainty in the relation be- 
tween SZ signal and mass is dominated by astrophysical pro- 
cesses and is only minimally dependent on cosmological pa- 
rameters given WMAP7 priors. 

We use the simulations discussed in Section [331 to deter- 
mine fiducial model values for A, B, and C. In particular, we 
use only the simulated maps containing the first-order thermal 
SZ component of the simulated clusters, without altering the 
maps by adding any noise or convolving with the ACT beam. 
Using the y value of the brightest 0.5' pixel for each cluster in 
the simulated thermal SZ maps, we solve for the best-fit val- 
ues of the three scaling relation parameters in Eq. [5] as well 
as the scatter in this relation, with a linear least squares fit. 
Figure [5] shows the simulated clusters along with the best-fit 
line, and Table [2] gives the corresponding best-fit values and 
errors of the fiducial scaling relation parameters in addition to 
the scatter S. The best-fit values for the B and C parameters 
are close to what we would expect from self-similar scaling 
relations as discussed in Appendix [Al 

There are a number of astrophysical mechanisms that could 
cause the observed relation between SZ signal and mass to 
differ from the fiducial relation. One is contamination of 
SZ decrements by radio or infrared galaxies at 148 GHz. 
Regarding radio galaxies, observations suggest that these 
gala xies show some pr e ference for residing in galaxy clus- 
ters (ICoble et al.1 120071: iLin & Mohrl l2^oTIOnetaiT 120091 



Cosmology from ACT Sunyaev-Zel'dovich Galaxy Clusters 



7 




3 5 10 15 20 30 

M 2m (lO l4 M e h- 1 ) 

FIG. 5. — Relation between true y and the cluster mass from simulations, 
including clusters with M 2 oo > 3 X 10 l4 M Q /i _I and z > 0.15. The best-fit 
scaling relation parameters in Eq.[5]are found with a least-squares fit, and 
the resulting best-fit surface is plotted in two-dimensions as the solid line 
above. The lognormal scatter between the true y values and the best-fit scaling 
relation is 26%. 

iMandelbaum et all l2009t) . However, using a model of ra- 
dio galaxies that describes their correlation with halos, the 
amount of contamination expected from rad io galaxies was 
found to be negligible for (Sehgal et al. 2010). For redshifts 
< 1, star formation, which is responsible for infrared galaxy 
emission, is expected to be quenched in high-density environ- 
ments. At low redshifts (z ~ 0.06) the fraction of al l galaxies 
that a re star forming galaxies is ~ 16% in clusters dBai et alJ 
2010). While this percentage is expected to increase at higher 
redshifts, given that the total infrared b ackground at 150 GHz 
is roughly 30 /iK (Fixse net al.lll998h . it is unlikely that in- 
frared galaxy contaminat ion could be sign ificant for clusters 
with y7cMB > 300 fiK. iLima et all (1201 Oh have also shown 
that the lensing of infrared galaxies by massive clusters should 
not introduce a significant bias in the measured SZ signals. 

Another way for the observed SZ signal to be lower than 
the fiducial model is if clusters have a significant amount of 
nonthermal pressure. This pressure would not be observed 
as part of the SZ signal, however, it would play an important 
role in counteracting the gravitational pressure from the clus- 
ter mass. Such nonthermal pressure can take the form of small 
scale turbulence, bulk flows, or cosmic rays. Simulations 
and observations suggest contributions to the total pressure 
from cosmic rays to be about 5- 10% (Jubelgas et all 120081 : 
iPfrom mer & EnBlin 2004 |) and from turbulent pressure to be 
between 5-20% (lLau etail 120091: iMeneghetti etail l20"lfi 
iBurns et ai]l2010l) . with only the latter work suggesting levels 
as high as 20% and that largely at the cluster outskirts. These 
processes have a much larger impact on lower mass clus- 
ters and groups where the gravitational potent ial is not strong 
enough to tightly bind the cluster gas (e.g., Battagli a et alJ 
120 lot IShaw et alJl2010t iTrac et alJl2010l) . However, for the 
massive systems considered here, this again is not expected to 
be a significant issue. One astrophysical process that can have 
a significant affect on the cluster y values is major mergers. 
We certainly have at least one in our sample (Bullet cluster), 
but note the extreme rarity of such objects in general. 



4.3. Cluster Likelihood Function 

In order to constrain cosmological parameters with our 
cluster sample, we construct a likelihood function specific for 
clusters, and we map out the posterior distribution to find 
marg i nalize d distributions for each parameter. We follow 
ICashl (119791) who derived the likelihood function in the case 
of Poisson statistics giving 

ln£ = lnP r ({«,}|{A,}) = J>/lnA,- A,). (6) 

!=1 

P r is the probability of measuring {«,-} given modeled counts 
{A,}. Here Nb is the total number of observed bins in SZ sig- 
nal - redshift space, and A, is the modeled number of clusters 
in the ith bin. We also take the bin sizes to be small enough 
so that no more than one observed cluster is in each bin. The 
modeled cluster count, A,, is a function of the SZ signal and 
redshift of the given bin (which we call y° hs and z obs ) as well 
as the set of cosmological parameters, {c,}. The modeled 
count is also a function of the parameters of the SZ signal - 
mass scaling relation (A,B,C,S) given in Eq.[3] since it is the 
abundance of clusters as a function of mass that is tied to cos- 
mology via the mass function . For th is work we use the mass 
function given in lTinker et al.1 (120081) . A derivation of the full 
cluster likelihood function used in this analysis can b e found 
in Appendix IB"1 This likelihood is given by Eq. IB8l and is a 
function of the parameters {cj} and A,B,C,S. 

We assume normal errors of 2.2 x 10~ 5 ony obs (correspond- 
ing to an error in 3>7cmb of 60 /iK) and 0. 1 on z obs . We take 0. 1 
as the redshift uncertainty for convenience even though six of 
our clusters have spectroscopic redshifts. However, the red- 
shift error does not dominate the uncertainty of our results. 
We also assume Gaussian priors on A,B,C, and S centered 
around the fiducial values given in Table [2] with conservative 
1(7 uncertainties of 35%, 20%, 50%, and 20% respectively of 
the fiducial values. These priors were determined by find- 
ing the relation between SZ signal and mass from simulated 
thermal SZ maps with varying gas models. In particular, we 
use two simulate d the rmal SZ maps analogous to those dis- 
cussed in Section l4~2l with the gas physics models in these 
maps base d on the ad iabatic and the nonthermal20 models de- 
scribed in lTrac et al.l (12.010b . The adiabatic model assumes no 
feedback, star-formation, or other nonthermal processes that 
could lower the SZ signal as a function of mass. The nonther- 
mal model assumes more star-formation than the fiducial 
model and 20% nonthermal pressure support for all clusters 
at all radii, which is a larger amount of nonthermal pressure 
than generally suggested by X-ray observations and hydro- 
dynamic simulations (e. g., lLau et alj|2009t fMeneghet ti et alJ 
2010: IBurns et alJl2010l) . These two models span the range 
of plausible gas models for massive clusters given current ob- 
servations, and the la priors on the scaling relation parame- 
ters given above are generous given the range in parameters 
spanned by these models. 

4.4. Parameter Constraints 

The likelihood function described above was made into a 
standalone code module which was then int erfaced with the 
Mark ov chain software package CosmoMC dLewis & Bridld 
120021) . Usi ng CosmoMC, we r un full chains for the WMAP7 
data alone dLarson et alj |2010) and for the WMAP7 data plus 
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FIG. 6. — Likelihood contour plots of w versus erg showing la and 2<r marginalized contours. Left: Blue contours are for W MAP7 alone, and red contours are 
for WMAP7 plus ACT SZ detected clusters, fixing the mass-observable relation to the fiducial relation given in Section l4~2l Right: C ontou rs are the same as in 
the left panel, except that the uncertainty in the mass-observable relation has been marginalized over within priors discussed in Section l4~3l 



TABLE 2 

Best-fit Scaling Relation Parameters 



Model and Data Set A B C S 



Simulation Fiducial Values (5.67 ± 0.05) X 10~ 5 1 .05 ± 0.03 1 .29 ± 0.05 0.26 
wCDM WMAP7 + ACT Clusters (8.77 ± 3.77) X 10" 5 1.75±0.28 0.97±0.68 0.27±0.13 



TABLE 3 

COSMOLOGICAL PARAMETER CONSTRAINTS FOR <T 8 AND W 



Model and Data Set erg w 



wCDM WMAP7+BAO+SN 0.802 ± 0.038 -0.98 ± 0.053 

wCDM WMAP7 0.835±0.139 -1.11 ±0.40 

wCDM WMAP7 + ACT Clusters (fiducial scaling relation) 0.82 1 ± 0.044 -1 .05 ± 0.20 

wCDM WMAP7 + ACT Clusters (marginalized over scaling relation) 0.851 ±0.115 -1 . 14 ± 0.35 



our ACT cluster subsample. We assume a wCDM cosmo- 
logical model which allows w to be a constant not equal to 
-1, assumes spatial flatness, and which has as free parame- 
ter s: £lbh 2 , £l c h 2 , 0*, t , w, n s , ln[10 I0 A s ], and Asz as defined 
in lLarson et al.l (12010I) . The parameter erg is derived from 
the first seven of these parameters and is kept untied to Asz 
as the link between the two is in part what we are investi- 
gating. We run the WMAP7 plus ACT clusters chain under 
two cases: one where the values of A,B, C, and S are fixed 
to the fiducial values given in Section l4~2l and listed in Table 
p2] and one where A,B,C, and S are allo wed to vary within 
the conservative priors given in Section 14.31 For the latter 
case, we add these four new parameters to the CosmoMC 
code. At each step of the chain, CosmoMC calls the soft- 
ware package CAMB0 to generate both the microwave back- 
ground power spectrum and matter power spectrum as a func- 
tion of the input cosmology, and then the natural logarithms 
of both the WMAP and cluster likelihoods are added. We 
determine the posterior probability density function through 
the Markov chain proc ess and use a simple R-l statistic 
dGelman & Rubinll 19921) of R-l < 0.01 to check for conver- 

6 www.camb.info 



gence of the chains. 

The best-fit marginalized la and 2er contours, obtained 
from this process, are shown in Figure [6] for w and <jg. 
The blue contours show the constraints for WMAP7 alone, 
while the red contours show the constraints from the union 
of WMAP7 plus our ACT cluster subsample. The left panel 
shows the best-fit contours with the SZ signal - mass scaling 
relation fixed to the fiducial relation obtained from the sim- 
ulations. The right panels show the constraints allowing the 
four parameters of the scaling relation to vary. Table[3]lists the 
best-fit parameter values for erg and w with their 1 a marginal- 
ized uncertainties. Table [2] lists the best-fit scaling relation 
values and la uncertainties as well as the fiduc ial v alues ob- 
tained from simulations as discussed in Section [4721 for com- 
parison. We note that for the remaining seven parameters fit 
in the analyses combining WMAP7 plus ACT clusters (ftbh 2 , 
Vl c h 2 , f?„, r, n s , ln[10 10 Aj], and Asz), we find best-fit values 
consistent with the best-fit values from WMAP7 alone with a 
modest improvement in the marginalized errors. 

4.5. Stacked SZ Signal 

We also perform a stacking analysis of the nine clusters 
listed in Table 1 to measure average cluster SZ profiles, which 



Cosmology from ACT Sunyaev-Zel'dovich Galaxy Clusters 



9 



-50 



^-100 



E— 

< 



-150 



-200 




_l I I l_ 



_L 



J I I l_ 



_L 



J I I l_ 



5 10 
Radius (arcmin) 



15 



FIG. 7. — Average profile from stacking the 9 clusters presented in Table 
[T] (black solid line) compared with the average profile stacking 40 clusters 
with unsmoothed >'7cmb values greater than 300 fiK in simulated thermal SZ 
maps convolved with the ACT beam (dashed blue line). The dashed blue line 
shows the average profile for clusters simulated with the fiducial SZ model, 
while the dotted red (bottom) and dot-dashed green (top) lines show the same 
assuming the adiabat ic an d nonthermal20 SZ models respectively which are 
discussed in Section |4.3| For the profiles of the 9 clusters in the data, we 
removed a mean background level from the profile of each cluster. Error bars 
for the simulated clusters have been offset by 0. 1 ' for clarity, and are smaller 
than those from the data. Error bars for the adiabatic and nonthermal20 mod- 
els are not shown, but are of similar size as for the fiducial model. 



can be compared with simulations. We stack the 9 clusters in 
the data map prior to any filtering, after subtracting a mean 
background level for each cluster profile using an annulus 15' 
from the center of each cluster and 0.5' wide. The stacked av- 
erage profile is given by the solid black line in Figure [7] The 
same procedure is preformed on all simulated clusters with 
unsmoothed yTcMB values greater than 300 /iK in simulated 
thermal SZ maps. There are 40 of these simulated clusters 
in total over 6 different 455 square degree maps spanning the 
same redshift range as the data. These simulated clusters are 
stacked in thermal SZ maps convolved with the ACT beam 
to mimic the data, and their average profile is given by the 
dashed blue line in Figure[7j Error bars represent the standard 
deviation of the mean in each radial bin. The blue dashed line 
represents the stacked profiles of simulated clusters assuming 
the fiducial SZ model. The red dotted and green dot-dashed 
lines show the stacked profiles of simulated clusters assum- 
ing the adia batic and nonthermal20 SZ models discussed in 
Section l4~3l The error bars have not been included for the lat- 
ter two models in Figure [7j but they are of similar size as for 
the fiducial model. We find good agreement in the average 
profiles of the clusters in the data and simulated with the fidu- 
cial model as shown in Figure |7J which suggests that there is 
no significant misestimate of the SZ signal for these massive 
systems. 

5. DISCUSSION 

From Table [3] we see overall agreement between as and w 
as measured with only WMAP7 and as measured with the 
high-significance ACT cluster sample plus WMAP7. We find 
o-g = 0.821 ± 0.044 and w = -1.05 ± 0.20 if we assume the 



fiducial scaling relation, a decrease in the uncertainties on 
these parameters by roughly a factor of three and two re- 
spectively as compared to WMAP7 alone. This indicates 
the potential statistical power associated with cluster mea- 
surements. Marginalizing over the uncertainty in this scaling 
relation, we find cr 8 = 0.851 ±0.115 and w = -1.14 ±0.35, 
an uncertainty comparable to that of WMAP7 alone. We 
also see consistency when comparing these constraints to the 
best-fit constraints from WMAP7 plus baryon acoustic oscil- 
lations plus type la supernovae, which give ag = 0.802 ±0.038 
and w = -0.980 ± 0.0 53 for a wCDM cosmological model 
( Kom atsu et alJ l2QToh . As the latter are all expansion rate 
probes, this suggests agreement between expansion rate and 
growth of structure measures. Both also show w is consistent 
with -1, giving further support to dark energy being an energy 
of the vacuum. 

These results are also consistent with analyses from X-ray 
cluster samples giving a s (fl m /0.25) 0A 1 = 0.813 ±0.013 (staf ) 
±0.024 (sys) and w = -1.14 ±0.21 ((yikhlinin et al. 2009), 
Q m = 0.23 ± 0.04, cr 8 = 0.82 ± 05, a nd w = -1.01 ± 0.20 
for a wCDM model (fMantz et alj|2010l) . and Q m = 0.30+g gl 
o-s = .85^ from WMAP5 plus X-ray clusters dHenrv et alJ 
2009) . We also find consistency with optical samples yield- 
ing o-«(fl m /0.25) a41 = 0.832 ±0.033 for a fl at ACDM model 
dRozo et al.l2010l) . IVanderlinde et alJd20Tob find cr 8 = 0.804 ± 
0.092 and w = -1.049 ± 0.291 for a wCDM model using SZ 
clusters detected by SPT plus WMAP7. 

This analysis also suggests consistency between the fiducial 
model of cluster astrophysics used here to describe massive 
clusters and the data. Table[3]shows agreement in best-fit cos- 
mological parameters between growth rate and expansion rate 
probes when we hold fixed our fiducial relation between SZ 
signal and mass. When we allow the scaling relation parame- 
ters to be free, we find best-fit values that are broadly consis- 
tent with those of our fiducial relation. We note that while the 
lcr range of the B parameter is higher than the fiducial value, 
the fiducial value is enclosed by the 2er range of 1 -75^q'1 The 
higher value of the B parameter may indicate some curvature 
in the true scaling relation away from the fiducial model at 
the high-mass end. This may also be suggested by Figure [5] 
where the simulated clusters seem to prefer higher y values 
than the fiducial relation would suggest for the most massive 
systems. The agreement between cosmological parameters 
from expansion rate and growth of structure probes when fix- 
ing the SZ signal - mass scaling relation to the fiducial model 
and the broad agreement between fiducial and best-fit scal- 
ing relation parameters when the latter are allowed to be free, 
suggest our data is broadly consistent with expectations for 
the SZ signal of massive clusters. This is also suggested by 
comparing the stacked SZ detected clusters in the data with 
simulations as shown in Figure|7j 

We would expect the above to be the case as massive clus- 
ters have been studied far better than lower mass clusters with 
a variety of multi-wavelength observations. In addition, a 
number of astrophysical processes that are not perfectly un- 
derstood, such as nonthermal processes and point source con- 
tamination, affect the gas physics of lower mass clusters much 
more than that of the most massive systems. In general, these 
processes tend to suppress the SZ power spectrum over that 
of a straightforward extrapolation based on the most massive 
systems. This is an important effect as lower mass systems 
(< 1O 14 M ) contribute as much to t he SZ power spectrum at 
/ ~ 3000 as systems at higher mass dKomatsu & Se liak 2002; 
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iTrac et ai1 l2010). The power spectrum n ear / = 3000 has bee n 
recently measur ed b y and discussed in | Lueker et al.l (120101) . 
iDas et all dloToh and lDunklev etafl d2010i) . 

There are a number of ways the cosmological constraints 
presented here could be further improved. Clearly the largest 
uncertainty is the relation between SZ signal and mass, and 
further X-ray observations of massive clusters, particularly at 
higher redshifts where X-ray observations have been limited, 
would help to calibrate this relation. Further targeted observa- 
tions of massive clusters at millimeter-wave frequencies with 
enough resolution and sensitivity to identify point sources 
would offer a better handle on contamination levels. In addi- 
tion, an analysis using multiple frequency bands, which would 
employ the spectral information of the SZ signal, may be 
helpful in determining cluster sizes and measuring integrated 
Ts. This could help reduce the scatter in the relation between 
SZ signal and mass. Spectroscopic redshifts of all the clusters 
in a given SZ sample would also help to reduce uncertainty 
on the cosmological parameters. In addition, millimeter-wave 
maps with lower instrument noise, would greatly reduce the 
scatter between the recover ed and true SZ signal . Such maps 
are expected with AC Tpol (Niemack et al. 2010) and SPTpol 
( McM ahon et al]l2009l) coming online in the near future. 

With continued SZ surveys such as ACT and SPT and 
their polarization counterparts, in addition to data forthcom- 
ing from the Planck satellite, we will no doubt increase the 
number of SZ cluster detections. We anticipate that upcom- 
ing larger galaxy cluster catalogs will make significant contri- 
butions to our understanding of both cluster astrophysics and 
cosmology. 
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APPENDIX 

A. SELF-SIMILAR SCALING RELATION BETWEEN SZ SIGNAL AND MASS 
If clusters were self-similar and isothermal, then we would expect the scaling relation between SZ signal and mass to be 

IU cx M 5 h ^E(z) 2 / 3 f g . ds /dl (Al) 

where Thaio is the Compton y-parameter integrated over the surface of the cluster in units of arcmin 2 , and E(z) = [^ m (l +z) 3 + 
J^a] 1 / 2 for a flat ACDM cosmology. The angular diameter distance is denoted by d&, and / gas is the gas mass fraction. For the 
Compton y-parameter integrated over a fixed aperture we have 

Taperture OC Thalo — , (A2) 

\ "halo / 

where /^aperture oc d& and T^aio oc M^J o E(z)~ 2 ^ 3 . Note that this equation is appropriate if the aperture size is smaller than the size 
of the cluster. The above gives 

^aperture OC M halo £(z) 2 / gas . (A3) 

To write Taperture as a function of (1 +z) we note that at z=0.5 (the mean redshift of our cluster sample) E(z) oc (1 + z ) - 835 for 
ft m = 0.27. Thus 

Aperture «M^ alo (l+z) C , (A4) 



where B= 1.0 andC= 1.67. 
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B. CLUSTER LIKELIHOOD FUNCTION 

Below we describe the construction of the likelihood function for SZ clusters detected in millimeter-wave surveys. From 
Poisson statistics, the probability of observing n,- counts expecting A,- counts is 

\ n, — \ t 

Pr(ni\\d = -!—-. (Bl) 

Given a data set of {«,-} counts in Nb observed bins and a corresponding prediction, {A,}, the probability of the data given the 
prediction is 

^(W|{A,}) = n^^, (B2) 

where A; = f > r (v obs ,z obs ,A,B,C,5,{c y })A^Ay obs Az obs . Here P r (y° bs ,z° hs ,A,B,C,S, {cj}) is the probability of observing a cluster in 
bin i, and N is a normalization factor giving A, units of counts (see below). The observed SZ signal and redshift of a given cluster 
are denoted by y obs andz obs , and Ay obs and Az obs denote the size of the bin. The parameters A, B,C, S describe the scaling relation 
between SZ signal and mass and are defined below. The cosmological parameters are indicated by {cj}. 
If we allow the bin sizes to be small enough that each observed bin holds no more than one observed cluster, then 

n N b 

\nPr({ni}\{\i}) = J2^i~J2 Xi = lnC 

i=l i=l 

as given in ICashl (119791) . Note that the ln(n,!) term has been dropped as it is independent of any change in parameters, and n 
represents the total number of clusters observed. Thus we have 

ln£ = S2ln(P r (y oh \z oh \A,B,C,S,{cj})Ndy° hs dz obs )- f f P r (y oh \z oh \A,B,C,S,{cj})Ndy oh *dz ohs (B4) 



with 



P r (y obs ,z obs ,A,B,C,S,{cj}) = JJJ W b ^^ ob ^A,B,C,5,{c ; },^^z ttu ^lnM ttue )^V^ true rflIlM tme (B5) 
P r (y obs , z° bs \A,B,C,S, {e,},/™ 6 , z trae , InM^PAA ,B,C,S, {cj},y tme , z trae , InM^dy^dz^dlnM 1 ™ (B6) 



P r (y ohs , z° bs |/ ue , z tme )P r (y tme \A,B,C,S, {cj} , z tme , \nM tme )P r {A,B ,C , S, { Cj }, z tme , \nM tme )dy l ™ dz l ™ d\nM aae (B7) 

/oo p oo poo 

dlnM true / dy tme / dz tme P r (y obs \y tme )P r (z ohs \z tme )P r (y tme \A,B,C,S,^ 
oo JO JQ 

x P r ({cj})P r (A)P r (B)P r (C)P r (S) (B8) 

using the definition of conditional probability. Here P r ({cj}) is any external prior on {cj} such as a WMAP prior. 
We assume the following SZ signal - mass scaling relation with log normal scatter, S, 



where M = 5x 1O 14 M /i 1 and zo = 0.5. This gives 

P fWA r r 9 ^« lnM true, 1 rTn / -(In/™ -fllnAf™ - On(l +z a n - lnA + B\nM + Cln(l +z )) 2 \ 

P r (y \A,B,C,S,z ,lnM )=^=__exp^ — 2 j. (BIO) 



We also assume Gaussian priors on the scaling relation parameters as indicated by simulations, giving 



1 f-(A-A ) 



^T^n^r'- <B11) 



Similar relations hold for P r (B), P r (C), and P r (S). 
From the mass function we have 
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v ' 11 1SJ d\aM Vme dz [me N 

where here n is the number density of clusters. N is the total number of clusters when the above mass function is integrated over 
d\nM tme andt/z^ 6 . 

We also assume for the uncertainty on the observed SZ signal and redshift that 

P r (y° hs \y lme ) = -7^— exp ^ / } (B13) 

1 / —( 7 °bs _ true\2 

_obs| true\ .„ t \& <■ ) 



P,.(z OBS |z true ) = -=- exp '~ - ' (B14) 
V2^cr z V 2cr z 2 / 

where these two expressions should also be multiplied by 1+erf(rlr „ 2 c ^ CT ) since the limits of integration are from to oo. 
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